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Abstract 

It is shown that under certain conditions the resonant transport in meso- 
scopic systems can be described by modified (quantum) rate equations, which 
resemble the optical Bloch equations with some additional terms. Detailed mi- 
croscopic derivation from the many-body Schrodinger equation is presented. 
Special attention is paid to the Coulomb blockade and quantum coherence 
effects in coupled quantum dot systems. The distinction between classical 
and quantum descriptions of resonant transport is clearly manifested in the 
modified rate equations. 
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I. INTRODUCTION 



Over the last decade a great interest has been paid to artificially fabricated nanostruc- 
tures containing discrete number of quantum states. The discreteness of quantum states 
manifests itself in peculiar transport properties of these systems as, for instance, in the 
Coulomb blockade oscillations [|T|. Actually, the study has been mostly concentrated on the 
quantum transport through single devices (quantum dots). In fact, more interesting quan- 
tum mechanical effects can be found in coupled nanostructures devices, where the quantum 
interference may strongly influence the resonance current. The impressive progress in micro- 
fabrication technology now allows to extend the experimental investigation to these systems 
too. For instance, the transport properties of coupled dots are presently under intensive 
study 

For description of quantum transport through single quantum dot (quantum well) the 
"classical" rate equations are usually used f^-^]. They can be derived either by using 
nonequilibrium Green's functions technique [0, or directly from the Schrodinger equation 
The situation is different for coupled wells with aligned levels. The quantum transport 
through these devices goes on via quantum superposition between the states in adjacent 
wells. It is thus quite obvious that non-diagonal density matrix elements would appear 
in the equations of motion. These terms have no classical counterparts, and therefore the 
classical rate equations have to be modifled. A plausible modiflcation of master equations 
for some particular cases of the resonance tunneling through double-dot structures has been 
proposed by Nazarov f^. A more general case, though without account of Coulomb inter- 
action, has been considered in |T^, where modifled rate equations have been proposed by 
using an analogy to the optical Bloch equations [|l^ . However, no microscopic derivation of 
the modifled rate equations has been presented yet. 

In this paper we derive the rate equations for a general case of resonant transport through 
mesoscopic systems, starting with the many-body Schrodinger equation, with special atten- 
tion being paid to the Coulomb blockade and coherent quantum mechanical effects. Our 



2 



main goals are, first, to substantiate and generalize the previously suggested rate equa- 
tions and second, to determine the region of validity of the rate equations for description of 
quantum transport in general. Also, we believe that the microscopic derivation of the rate 
equations will provide a better understanding of the correspondence between quantum and 
classical description of carrier transport in mesoscopic systems. 

The plan of the paper is the following. In Sect. 2 we give a detailed derivation of the 
transport rate equations through a single quantum well (dot). In order to present our 
method most lucidly, we neglect in this section the Coulomb interaction and spin effects. 
These effects are considered in Sect. 3. In Sect. 4 we derive the modified rate equations for 
coupled well structures, taking into account the Coulomb and spin effects. An example of 
coherent resonant transport with inelastic transitions is studied in Sect. 5. The derivation of 
rate equations performed in this case allows to establish their correct form valid in general 
case of quantum transport. The general case and an example of coherent resonant transport 
with inelastic transitions in the presence of strong Coulomb blockade are presented in Sect. 6. 
The last section is a summary. 

II. SINGLE- WELL STRUCTURE 

Let us consider a mesoscopic "device" consisting of a quantum well (dot), coupled to two 
separate electron reservoirs. The density of states in the reservoirs is very high (continuum). 
The dot, however, contains only isolated levels. We first demonstrate how to achieve the 
reduction of many-body Schroedinger equation to the rate equation in the simplest example. 
Fig. 1, with only one level, Ei, inside the dot. We also ignore the Coulomb electron-electron 
interaction inside the well and the spin degrees of freedom. Hence, only one electron may 
occupy the well. With the stand simplifications, the tunneling Hamiltonian of the entire 
system in the occupation number representation is 

?i = ^ Eia\ai + Eittitti + ^ E^alar + ^ r2;(a|ai -|- a\ai) + ^ r2^(ajai -|- a|a^) . (2.1) 

I r I r 
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Here the subscripts / and r enumerate correspondingly the (very dense) levels in the left 
(emitter) and right (collector) reservoirs. For simplicity, we restrict ourselves to the zero 
temperature case. All the levels in the emitter and the collector are initially filled with 
electrons up to the Fermi energy Ep and Ep, respectively. This situation will be treated as 
the "vacuum" state |0). 

This vacuum state is unstable; the Hamiltonian Eq. (p.l| ) requires it to decay exponen- 
tially to a continuum state having the form a|a;|0) with an electron in the level Ei and 
a hole in the emitter continuum. These continuum states are also unstable and decay to 
states aja; |0) having a particle in the collector continuum as well as a hole in the emitter 
continuum, and no electron in the level Ei. The latter, in turn, are decaying into the states 
a|aja;a;/ 10) and so on. The evolution of the whole system is described by the many-particle 
wave function, which is represented as 



bo{t) + '^bii{t)a\ai + '^bir{t)alai + ^ biii'r{t)a\alaiai' + . . . 

I l,r l<l',r 



|0), (2.2) 



where b{t) are the time-dependent probability amplitudes to find the system in the corre- 
sponding states described above with the initial condition 6o(0) = 1, and all the other 6(0)'s 
being zeros. Substituting Eq. ( ^^ into in the Schrodinger equation i|\I/(t)) = 7i|\I'(t)), 
results in an infinite set of coupled linear differential equations for the amplitudes b{t). 
Applying the Laplace transform 



oo 

iEti 



b{E) = / e'''%{t)dt (2.3) 
Jo 

and taking account of the initial conditions, we transform the linear differential equations 
for b{t) into an infinite set of algebraic equations for the amplitudes b{E), 

Eb,{E) -Y,nibu{E) = I (2.4a) 

{E + Ei- Ei)bii{E) - QME) -J2Qrbir{E) = (2.4b) 

r 

{E + Ei- Er)biriE) - QMe) -Y.^Au'r{E) = (2.4c) 

I' 

{E + Ei + El, -El- Er)bm'r{E) - niblr{E) + Vtibvr{E) - ^r'bll'rr'{E) = (2.4d) 
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Eqs. (|2.4| ) can be substantially simplified. Let us replace the amplitude b in the term 
J2 of each of the equations by its expression obtained from the subsequent equation. For 
example, substitute bii{E) from Eq. (|2.4b| ) into Eq. ( p.4a| ). We obtain 

e-y: — 

I E + El — El 

Since the states in the reservoirs are very dense (continuum), one can replace the sums over / 
and r by integrals, for instance J2i ^ I Pl{Ei) dEi , where Pl{Ei) is the density of states in 
the emitter. Then the first sum in Eq. ( p.5|) becomes an integral which can be split into a sum 
of the singular and principal value parts. The singular part yields —iQ{Ep + E ~ Ei) 
where = 2hpl{Ei)\Ql{Ei)\^ is the level Ei partial width due to coupling to the emitter. 
Let us assume that Ep ^ Ei ^ Ep, i.e. the bias is large and the energy level is deeply 
inside the band. In this case the integration over -variables can be extended to ±cx3. 
As a result, the theta-function can be replaced by one, and the principal part is merely 
included into redefinition of the energy Ei. Also, the second sum (integral) in Eq. ( p.5| ) 
proves to be negligibly small. Indeed, let us replace bir b{Ei, Er, E), and assume weak 
energy dependence of Q on -Ez(r). Then one finds from Eqs. ( p.4| ) that the poles of the 
integrand in the i?/(i?r )-variable are on one side of the integration contour, and therefore 
this term vanishes. 

Applying analogous considerations to the other equations of the system (|2.4|) , we finally 
arrive to the following set of equations: 

{E + iTL/2)bo{E) =i (2.6a) 

{E + E1-E1+ iTR/2)bn{E) - nME) = (2.6b) 

{E + Ei~ Er + iVL/2)bir{E) - nrbu{E) = (2.6c) 

{E + Ei + Ev -Ei-Er + iTR/2)Kvr{E) - Qi~bir{E) + Qihr{E) = (2.6d) 



where Tr = 27ipR{Ei)\flR{Ei)\'^ is the level Ei partial width due to coupling to the collector. 



l.r 



E + El — El 



bir{E) 



I. 



(2.5) 
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Now we introduce the density matrix of the "device" . The Fock space of the quantum 
well consists of only two possible states, namely: \a) - the level Ei is empty, and \b) - 
the level Ei is occupied. In this basis, the diagonal elements of the density matrix of the 
"device", aaa and abb, give the probabilities of the resonant level being empty or occupied, 
respectively. In our notation, these probabilities are represented as follows: 

(Taa = Mt)\^ + T.Kit)\^ + E \bll'rr'{t)\^ + . . . 
l,r l<l',r<r' 

- + ^ii^ + + • • • , (2.7a) 

I l<l',r l<l'<l",r<r' 

^ ^i? + 41^ + + ■■■, (2.7b) 

where the index n in a^""^ denotes the number of electrons in the collector. The current 
/(t) flowing through the system is I(t) = eNji{t), where Nji{t) is the number of electrons 
accumulated in the collector, i.e. 

NR{t)=Y.n[aj-\t) + a!-\t)] (2.8) 

n 

The density submatrix elements are directly related to the amplitudes b{E) through the 
inverse Laplace transform 

^(n)^,) = E / ^~&^...■■■(i^)^^......(i^')e^(^'-^^* (2.9) 

l...,r... 

By means of this equation one can transform Eqs. ( |2.6| ) for the amplitudes b{E) into 
differential equations directly for the probabilities cr*^"^(t). Consider, for instance, the term 
(^il\t) = El Ibiiit)]"^, Eq. ( p7b| ). Multiplying Eq. (plGEQ by bli{E') and then subtracting the 
complex conjugated equation with the interchange E ^ E' we obtain 

/ ^(^' - ^ - '^^^ E hdE)bUEy^^^'-^^^ 

-J ^^2lmY^nMEK{Ey^'''-''^' = (2.10) 

One can easily deduce from Eq. ( [2. 91 ) that the first integral in Eq. ( |2.1U| ) equals to —i[&ll\t) + 
^R'^bbi^)]- Next, substituting 



^''^""^ = E + E,-E,+rTn/2 ^^-ll) 
from Eq. ( p.6b|) into the second integral of Eq. ( |2.10|) , and replacing X)i ^ / PL{Ei)dEi one 
can can perform the i?;-integration in the integral, thus obtaining iTiaj^Jit). Ultimately, 
Eq. ( p30|) reads a^\t) = TLCT^^J{t) - TRail\t). We can go on with this algebra for all the 
other amplitudes b{t). For instance, the above procedure applied to Eq. ( p.6d| ) converts it into 
a differential equation for the density- matrix element cr^^'', Eq. ( |2.7bD . The only difference 



with the previous example is an appearance of the "cross terms", like J2^ibi'r{E)Qiibi^{E'). 
Yet, these terms vanish after the integration over Ei(^r), just as the second term in Eq. (|2.5|) ). 
The rest of the algebra remains the same, so one obtains CTil\t) = TLa^^it) — r/jcr^^^(t). 
Finally we arrive to the following infinite system of the chain equations for the diagonal 
elements, and (Jbb\ of the density matrix, 

^!^it) = -TLcri^J{t), (2.12a) 

^^\t) = TLcrj^J{t)~Tnail\t), (2.12b) 

^^:J{t) = -TLcr^:Xt)+Tn4l\t), (2.12c) 

^il\t) = TLa^:J{t)-Tncril\t), (2.12d) 



Summing up these equations, one easily obtains differential equations for the total proba- 
bilities aaa = En c^ia^ and abb = En (^bb^- 

Oaa = -^L<^aa + ^R<ybb , (2.13a) 

abb = ^L(^aa - ^R(^bb , (2.13b) 

which should be supplemented with the initial conditions 

aaaiO) = 1, a,,(0) = 0. (2.14) 
Using Eqs. ( |2.8| ), ( |2.12| ) we obtain the total current 

I{t) = eNnit) = erR[aS\t) + all\t) + a^\t) + ...] = eTnabbit). (2.15) 



Thus the current I{t) is directly proportional to the charge density in the well. Solving 
Eqs. ( p. 131 ) and substituting ahtit) into Eq. ( |2.15| ), we obtain (for t oo) the standard 



formula for the dc resonant current, 



I/e=-^. (2.16) 

i- L + i- R 

Notice that whereas the time-behavior of the current I{t) depends on the initial condition, 
the stationary current / = I{t oo), Eq. (|2.16| ), does not. 

Equations ( p. 131) , derived from the many-body Schrodinger equation, coincide with the 



classical rate equations in the sequential picture for the resonant tunneling, obtained using 
nonequilibrium quantum statistical mechanics technique 0. In contrast, our approach starts 
directly from the many-body Schrodinger equation and will be straightforwardly extended 
to more complicated situations. Note, however, that the method can be applied only when 
the resonance energy is inside the band, and Tl,r <^ Ep — Ep. If the resonance is near 
the edge of the band, but the width of the resonance is much smaller than the band width, 
our method still can be applied, but only to the stationary case (t — > oo). Yet, the time- 
dependent Scrodinger equation cannot be reduced to the rate equations ( |2.13| ), and therefore 
this case is not a subject of this paper. 

III. COULOMB BLOCKADE 

Now we extend the approach of Sect. 2 to include the effects of Coulomb interaction. 
Consider again the quantum well in Fig. 1, taking into account the spin degrees of freedom 
(s). In this case the tunneling Hamiltonian (|2.1| ) becomes 



l,s s r,s 

+ J2^i{(4s(^is + aL«/s) + Y^r{alsais + aLors) + Ua\^aisa\_^ai_s , (3.1) 

l,s r,s 

where s = ±1/2, and U is the Coulomb repulsion energy. 
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Writing down the many-body wave function, in the occupation number rep- 

resentation, just as in Eq. (p.2|), and then substituting it into the Schrodinger equation 
z|\E'(t)) = 7i|\&(t)), we find a system of coupled equations for the amphtudes b{t) 

(3.2a) 
(3.2b) 
(3.2c) 



EboiE)-J2ni [b^-iiE) + b^iiE)\ =1 
I 

{E + Ei~ ESviE) - QME) - J2 ^I'bnii'iE) - ^ nXiE) = 

I' r 

{E + Ei- Er)bir{E) - n^E) -J2^i' [h'iE) + = 

I' 

{E + Ei + El, - 2Ei - U)bnii'{E) - Qib^iiE) - QibniE) 



(3.2d) 



In order to shorten notations we ehminated the index (1) of the level Ei in the amplitudes b, 
so that 6|(|)...(^) denotes the probability amplitude to find one electron inside the well with 
spin up (down), and the amplitude b^i,,,(t) is the probability amplitude to find two electrons 
inside the well. 

Eqs. ( |3.2| ) can be simplified by using the same procedure as described in the previous 
section. For instance, by substituting bir from Eq. ( |3.2c| ) and b^ui from Eq. ( |3.2d| ) into 
Eq. ( p.2b|) , and neglecting the "cross terms" on the grounds of the same arguments as in 
the analysis of Eq. (|2.5| )), we obtain 

^ PL{Ev)^\Ei.)dEi. pR{Er)n\Er)dEr 



E + El — E 



b^i{E)=Q (3.3) 



-oo E + Ei + Ev -2Ei-U JeR E + Ei - Er 

Since Ei Ei, the singular parts of the integrals in ( p.3|) are respectively —iQ{Ep + E — 
El + U) V'J2 and -iQ{E + Ei - E§) Tn/2, where 

TLiR) = 27TpL^n){Ei)\nLiR){Ei)\', r^(«) = 27rpi(ii)(^i + U)\nLiR){Ei + U)\\ (3.4) 

Here Pl{r) is the spin up or spin down density of states in the emitter (collector), Pl{r) = 
Pl{r)] = Pl{r)i- As in the previous section, we assume the resonance level being deeply 
inside the band, Ep <ti Ei Ep. If, in addition, Ei + U <^ Ep, the theta-function in 
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the singular parts of the integrals in (|3.3| ) can be replaced by one. In the opposite case, 
-El + f/ ^ Ep, the corresponding singular part is zero. 



Proceeding this way with the other equations of the system (|3.2|), we finally obtain 



{E + iTL)bo{E) =1 (3.5a) 

(E + Ei-E, + iT'j2 + iTr/2%i{E) - nME) = (3.5b) 

{E + Ei-Er + iTL)bir{E) - QrME) = (3.5c) 

{E + Ei + Ev -2Ei-U + 21^6^0' {E) - ^ih^v {E) + nih^i {E) = (3.5d) 



Eqs. (oTo) can be transformed into equations for the density matrix of the "device" by 



using the method of the previous section. Since the algebra remains essentially the same, 



we give only the final equations for the diagonal density matrix elements cr^^t), al]^}{t), 



cr^^l(t) and cr^"-'(t). These are the probabilities to find: a) no electrons inside the well; 
b) one electron with spin up (down) inside the well, and c) two electrons inside the well, 
respectively. The index n denotes the number of electrons accumulated in the collector. We 
obtain 

^i? = -2r,aW + Tnal-'^ + T^a^f (3.6a) 

<^ = -(r^ + r^)a(:) + r^aW + r^^a^r^) (3.6b) 

<^ = -(r^ + rn)ai:l + r^ai^) + r^a^r^) (3.6c) 



These rate equations look as a generalization of the rate equations ( |2.12| ), if one allows the 
well to be occupied by two electrons. The Coulomb repulsion leads merely to a modification 
of the corresponding rates F —>■ T', due to increase of the two electron energy. 

Summing up the partial probabilities we obtain for the total probabilities, a{t) = 
Y.n (^^""Ki) ^ the following equations: 

CTaa = -'^TiCTaa + ^R'^bb't + ^R(^bbl (3.7a) 
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0"b6T = "(^L + rij)cr5f,| + TiCTaa + ^'r<^cc (3.7b) 

abbi = -(r^ + r/j)afeb| + TLaaa + r'^o-cc (3.7c) 

&CC = -2r'^crcc + r'^o-ftbi + T'^dhbi, (3.7d) 

and for the current 

I(t) = ^n[(T(")(t)] = er^ [atbrit) + abbM + 2eT'j,a,,it) (3.8) 

n 

Eqs. ( p.7|) , can be solved most easily for dc current, / = I{t oo). In this case 
cr = 0, and Eqs. (|3.7|) turn into the system of linear algebraic equations. One also finds 
from Eq. ( p.7|) that aaa + crbb^; + cTfefei + cTcc = 1 • The latter implies that dc current does not 
depend on the initial conditions. Finally we obtain 

i/e = 2r,r^(r- + r^) 

If El < < El + f/, one finds from Eq. (|0|) that = 0, so that the state with two 
electrons inside the well is not available. In this case one obtains from Eq. ( p.9| ) for the dc 
current 

I/e=:p^ (3.10) 



It is interesting to note that this result is different from Eq. ( 2.16| ), although in both cases 



only one electron can occupy the well. However, if the Coulomb repulsion effect is small, 
i.e. = Tl ji, Eq. (|3.9| ) does produce the same result as Eq. (|2.16| ), provided the density 
of states is doubled due to the spin degrees of freedom. 

One can also consider the case when the Fermi level in the right reservoir Ep lies above 
the resonance level Ei, but below Ei + U, so that Tji = 0, Eq. (p.3|). Then the resonant 
transitions of electrons from the left to the right reservoirs can go only through the state 
with two electrons inside the well. Using Eq. (|3.9|) one finds for the dc current 



which coincides with the result found by Glazmann and Matveev 
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IV. DOUBLE- WELL STRUCTURE 



A. Non-interacting electrons. 



Now we turn to the coherent case of resonant tunnehng. Let us consider the coupled- 
well structure, shown in Fig. 2. We assume that both levels E12 are inside the band, i.e. 
Ep <^ El, E2 Ep. In order to make our derivation as clear as possible, we begin with the 
case of no spin degrees of freedom and no Coulomb interaction. The tunneling Hamiltonian 
for this system is 



= ^ Eia\ai + Eia\ai + i?202'22 + X! ^rCilar 

+VLQ{a\a2 + a^ai) + ^ Vti{a\ai + a\ai) + ^ ^^(0^,02 + 4 



2"'r ) 



(4.1) 



where ajg, 0,1,2 are creation and annihilation operators for an electron in the first or the 
second well, respectively. All the other notations are taken from Sect. 2. The many-body 
wave function describing this system can be written in the occupation number representation 



as 



Lr 



I W l<l',r 



|0), 



(4.2) 



Substituting Eq. ( ^^ ) into the Shrodinger equation with the Hamiltonian ( |4.1| ) and per- 
forming the Laplace transform, we obtain an infinite set of the coupled equations for the 
amplitudes b(t): 

Ebo{E)-J2nibu{E) = t (4.3a) 
{E + Ei- Ei)bu{E) - nME) - nob2i{E) = (4.3b) 
{E + Ei- E2)b2i{E) - nobu{E) - J2 ^I'bmi'iE) - J2 ^ME) = (4.3c) 

I' r 

{E + Ei + El, -El- E2)bi2w{E) - Qib2i{E) + nib2i'{E) -Y^nMvr{E) = (4.3d) 
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Using exactly the same procedure as in the previous sections, Eq. (|2.5| ), we transform Eqs. 
( [4 .31) into the following set of equations: 

{E + iTl/2ME) = I (4.4a) 

{E + Ei- Ei)bu{E) - nME) - nJ)2i{E) = (4.4b) 

{E + E1-E2 + iTl/2 + iTR/2)h2i{E) - QobiiiE) = (4.4c) 

(E + Ei + Ev -E,-E2 + tTn/2)buii'iE) - ^vhiiE) + ^ihriE) = (4.4d) 



The amplitudes h{t) determine the density submatrix of the system, cr^"\ in the corre- 
sponding Fock space: \a) - the levels £'1^2 are empty, \h) - the level Ei is occupied, |c) - the 
level E2 is occupied, \d) - the both level Ei^2 are occupied; the index n denotes the number 
of electrons in the collector. The matrix elements of the density matrix of the "device" can 
be written as 



0"aa 


/ y aa 
n 


^\bo{t)\' + T.K{t)\'+ E \bii'rAt)\' + . . . 

l,r l<l' ,r<r' 




(4.5a) 




\ " (n.) 

n 


= E l^lK^)!' + E \blll'rit)\' + E IblU'l'Wr 
I l<V,r l<l'<l",r<r' 


{t)n... 


(4.5b) 




n 


= EI^2K^)P+ E l^2«/'r(t)P+ E \b2Wl"rr 
I l<l',r l<l'<l",r<r' 


(t)p+... 


(4.5c) 




\ ^ (n) 

n 


= E 1 ^12«' it)\^+ E 1 bl2U'l"l"'rr' (t) T + • 
Kl' l<l'<l"<l"',r<r' 




(4.5d) 




= E^6c 

n 


^T^biimiit) + E bm'r{t)b*2u>rit) + ■■■ 

I l<l',r 




(4.5e) 



Now we transform Eqs. ([4.4| ) into differential equations for a^^\t). Consider for in- 
stance the term all^ = J2i \bii{t)\'^ , Eq. ( [4.5b| ), where the amplitudes bu are determined 
by Eq. ( |4.4b| ). Multiplying Eq. ( |4.4b| ) by b\i{E') and subtracting the complex conjugate 



equation with E ^ E', we find 
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J:{E' - E)h{E)bl,{E') bl{E')bu{E) - HE)h\,{E')] 

I I 



After applying the inverse Laplace transform, Eq. 
becomes —i&^^\t). Next, substituting 

~ nME) + nohi{E) 



- noJ2r2iiE%iiE) - UE)bl,iE')] = (4.6) 
the first term in this equation 

(4.7) 



E + Ei-E^ 

from Eq. ( [4.4b| ) into the second term of Eq. ( |4.6| ), and replacing the sum by an integral 
over El, we reduce this term to iV ]J)Q{E)bQ{E') . After the inverse Laplace transform it 
becomes iTia^J^t). Notice that the "cross term", oc ^7o^^^o^2/! does not contribute to the 
integral over Ei, since the poles of the integrand in the i^^-variable lie on one side of the 
integration contour (cf. the second term of Eq. (|2.5|) ). The third term of Eq. ( [4.6|) turns to 
be no[ai''J{t)-aS,\t)], after the inverse Laplace transform. Finally we obtain a differential 



equation for the density submatrix element all\ 

4?W = r.ai^+^f^o(ai?-aS) 



(4.8) 



In contrast to the rate equations of the previous sections, the diagonal matrix element abb 
is coupled with the off-diagonal density matrix element abc- 

The corresponding differential equation for abc can be easily obtained by multiplying 
Eq. ([4.4b|) by (-^') with subsequent subtracting the complex conjugated Eq. (|4.4c| ), mul- 
tiplied by bii. Afterwords by integrating over Ei we obtain 



Eventually we arrive to the following set of equations for a^"^^ 



(0) 



.(0) 
be 



a 



(n 



a 



a 



(n 



a 



a. 



(n-l) 



be 



a 



cb J ) 



-rrf-r,a(:?)-.fio(4?-a^:)) 

^{E2 - E,)al:^ + znoial-^ - ai:^] 
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^{rL + rR)a, 



(n) 
be 



(4.9) 

(4.10a) 
(4.10b) 
(4.10c) 
(4.10d) 
(4.10e) 



Using Eqs. ( |4.10|) we can find the charge accumulated in the collector, Nuit), and subse- 
quently, the total current, eN(t), as given by 

J(t)/e = N{t) =J2n + + cri:\t) + aSit)] = r« [a^t] + cr,M (4.11) 

n 

As in the previous examples, the current is proportional to the total probability of finding an 
electron in the well adjacent to the right reservoir. The off-diagonal elements of the density 
matrix do not appear in Eq. (|4.11| ). 

Summing up over n in Eqs. ( [4.10|) , we obtain the system of differential equations for the 
density matrix elements of the device 

C^aa = -'^L'^aa + ^R(^cc , (4.12a) 
abb = ^L<yaa + ^R<ydd + i^o{<^bc " cb) , (4.12b) 
C^cc = -^R<ycc - ^L(^cc - i^o{<ybc - CTcb) , (4.12c) 

CTdd = -r/jCTrfrf + r^cTcc , (4.12d) 
&bc = i{E2 - Ei)crbc + i^o{(Tbb - (Jcc) - ^(E^ + TR)abc , (4.12e) 

Eqs. ( [4.12| )) resemble the optical Bloch equations [|11|]. Note that the coupling with the 
reservoirs produces purely negative contribution into the non-diagonal matrix element's 
dynamic equation, Eq. ( [4.12e| ), thus causing damping of this matrix element. 

Eqs. (|4.12|) are solved most easily for the stationary current, / = I{t —>■ oo). Using 



(^aa + CTbb + CTcc + CTrfd = 1, WC obtaiu 

I/e=( ^^^^ ^ ^ (413) 

where e = E2 — E1. This result coincides with the one found in the framework of one-electron 
approach [|1^,0]. 



B. Coulomb blockade. 



The extension of the rate equations ( ^.12| ) for the case of spin and Coulomb interaction 
is done exactly in the same way as in Sect 3. Here also the rate equations for the device 
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density matrix are obtained only for £"12 + U being inside or outside the band, but not close 
to the band edges [Ep -C -Ei,2 + U <^ Ep ot Ei^2 + > Ep). Eventually we arrive to the 
rate equations of type Eqs. ( |4.12| ), but with the number of the available states of the device 
changed due to additional (spin) degrees of freedom and Coulomb blockade restrictions. The 
Coulomb repulsion manifests itself also in a modification of the transition amplitude Q and 
the rates T's, Eq. (gj). 

In the case of large Coulomb repulsion, some of electron states of the device are outside 
the band (the Coulomb blockade). As a result, the number of the equations is reduced. 
Consider, for instance, the situation where the Coulomb interaction U of two electrons in 
the same well so large that E12 + t/ ^ Ep, but the Coulomb repulsion of two electrons in 
different wells, U, is much smaller, so that Ei^2 + U Ep. Then the state of two electrons 
in the same well is not available, but two electrons can occupy different wells. In this case 
the rate equations for the corresponding density matrix elements of the device are 

aaa = -2TL(Taa + ^R{<ycc\ + <ycci) , (4.14a) 
0"66T = ^L<yaa - r^((Trfrf-f| + a ddu) + ^^o(o"bcT " Cy cb]) , (4.14b) 

c^ccT = -r/jCTcct - 2r^(TccT - "^^0(0"^ - o-c6t) ) (4.14c) 
c^ddTT = -r'/jCTddi^ + r'^CTccT , (4.14d) 
dfecT = i{E2 - Ei)abc^ + if^o(o-66T " ^cct) - ^C^^'l + ^R)(^bcr , (4.14e) 

where r^^-^^ = 2np^ji){Ei + U)\QL(^ji){Ei + f/)p. Here for the shortness we wrote only 
the equations for the "spin up" component of the density matrix. The same equations are 
obtained for the "spin down" components of the density matrix. The total current is 

I/e = Tn^acc] + cTcci) + ^'nif^ddn + (^ddu + ^ddu + (^ddu)- (4.15) 

It is quite clear that the "spin up" and "spin down" components of the density matrix 
are equal, i.e. abb] = o'bbi = o'bb, the same holding for dec, Cdd components. Therefore Eqs. 
( [4.14| ), ( [4.15[ ) can be rewritten as 
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Gaa = -2TLO-aa + '^^RCTcc , (4.16a) 

&bb = ^LCTaa + 21^(7^^ + iQo{abc - (^cb) , (4.16b) 

O-cc = -'^RO-cc - 2r'^ (Tec - i^oicTbc - CTcb) , (4.16c) 

CTdd = -r^CTrfrf + r'^fTcc , (4.16d) 

abc = i{E2 - Ei)abc + i^o{crbb - c^cc) - \ i'^^'L + ri?)crfec , (4.16e) 



and 



I/e = 2TRa^ + ArR(jM (4.17) 
Using aaa + 2crfe6 + 2crcc + 4(7^^ = 1 we obtain for the dc current 

° (2r^ + r«)2 ^ 2 ^ (2r^ + r^)2 

where e = E2 — Ei. Notice that the current ( [4.18|) differs from that given by Eq. ( [4.13| ) even 



for r'j^ = Tl and = F/j, despite the fact that in the both cases only one electron can 
occupy each of the wells. 



It is interesting to compare our result with that of Stoof and Nazarov [|I^ for the case of 
strong Coulomb repulsion between two electrons in different wells (-£^1,2 + U ^ Ep), where 
only one electron can be found inside the system. It corresponds to F^ = 0. In this case the 
dc current given by Eq. ([4.1 8|) is 



This result is slightly different from that obtained by Stoof and Nazarov (by the factor 
two in front of F^). The difference stems from the account of spin components in the rate 



equations, which has not been done in WA 



V. INELASTIC PROCESSES 



As an example of a system with coherent tunneling accompanied by inelastic scattering, 
let us consider the coupled-dot structure shown in Fig. 3. In this system a resonant current 
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flows due to inelastic transition from the upper to the lower level in the left well. For sim- 
plicity, we restrict ourselves to non-interacting spin-less electrons. The Coulomb interaction 
and the spin effects can be accounted for precisely in the same way as we did in the pre- 
vious sections, namely by allowing for states with doubly occupied levels (excluding states 
violating Coulomb restrictions) and modifying transition amplitudes and inelastic rates. 
The tunneling Hamiltonian of the system has the following structure 

'H = ^ Eia\ai + Eia\ai -\- ii^2a2«2 + E^alas + ^ E^^cl^Ca + ^ Erolar + l^o(a2«3 + «3«2) 

/ a r 

+ ^Qi{alai + a{ai) + ^Q^{alaicl^ + a\a2Ca) + ^Qr{al(^3 + 4«r)- (5.1) 

I a r 

Here the subscript a enumerates the states in the phonon bath and fl^ is the corresponding 
coupling. The many particle time-dependent wave function of the system is 



mt)) = 



^o(^) + J2 bii{t)a\ai + b2ia{t)alaicl + ^ b3ia{t)alaicl+ 

I l,cx 



+ h2U'a{t)a\alaiai>cl+ ^ bi3ii'a{t)a[alaiai>cl + . . . 

l<l',a l<l',a 



|0). (5.2) 



Repeating the procedure of the previous sections we flnd the following set of equations for 
the Laplace transformed amplitudes, b{E): 

(E + irL/2)bo^i (5.3a) 

{E + E1-E1 + iTj2)bn - ^ibo = (5.3b) 

(E + Ei-Ea-E2 + iVL/2)b2ia - K^'bii - ^ohia = (5.3c) 

{E + Ei-Ea-Es + irL/2 + irR/2)b3ia - ^ohia = (5.3d) 

{E + El + Ell — El — E2 — Ea)bi2ll'a — ^l'b2la + ^lb2l'a — ^obisii'a — (5.3e) 

{E + Ei + El, -Ei-E3-Ea + iTj2 + iTR/2)bmi'a 

— ^obl2ll'a — ^I'bsia + ^ibsi'a — (5.3f) 



where Tin = '^'^'Pphl^^'^l'^ is the partial width of the level Ei due to phonon emission and Pph 
is the density of phonon states. 
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The density matrix elements of the device is (Jij{t) = J2n'^ij\t)^ where alj\t), are 
related to the amplitudes b{E) via Eq. ( p.9|) . All possible states electron states of the device 
are shown in Fig. 4. Using the previous section procedure for diagonal matrix elements 



we obtain master equations analogous to Eq. ( 4.12 ), in which transitions between isolated 



levels E2 and E3 take place through the coupling with non-diagonal matrix elements. These 



equations have the appearance of the optical Bloch equation |rT|. However, the master 
equation for the non- diagonal matrix element, o"e/, contains an additional term. Therefore, 
we present the derivation of the master equations for "coherences" cTg/ and acd in some detail. 
Consider for example the non-diagonal density submatrix elements a^J = J2i,ab2ia{t)blia{t) 
and (Tgj = J2i<i'.abi2ii'a{t)bl^ii,^(t). The differential equation for cr^°''(t) can be obtained by 
multiplying Eq. ( ^.3c|) by b^i^i^E') with subsequent subtraction of the complex conjugated 
Eq. ( |5.3d| ), multiplied by h2ia{E). Then using Eq. ( p.9|) , we obtain 

^1? = ^(^3 - E,)af] + ^fio(a(? - aS) - h2T, + r^)^^? . (5.4) 



Similarly, multiplying Eq. ( p.3e| ) by h\^ii,^{E') and Eq. ( p.31| ) by bi2ii'a{E), we find the dif- 

(0) 
ef 

af) = i{E, - E2)af) + zn,{a^^^ - afj) - ^(r„, + Tn)a^^j/ - ^A. (5.5) 



ferential equation for a^) (t) 



where 



-h,2Wa{E)nvhli^{E') + h,2n'a{E)nihli,^{E')\ e^^^'"^)* (5.6) 



Substituting the amplitudes buwa from Eq. ( |5.3e|) and bl^n,^ from Eq. (|5.3i| ) into Eq. {f).6\), 
and replacing the sum over /(/') by the corresponding integral, we find —iA = r^cr^ll^ It 
implies that the non-diagonal density matrix a^j given by Eq. ( |5.5| ), is coupled with acd via 
a single electron transition from the emitter to the left well. Such a term does not appear 
in the Bloch equations, which deal with two-level systems. 
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Summing up over n in the rate equations for the density submatrix (t) we obtain the 

set of rate equations for the density matrix of the device 

O-aa = -'^LO'aa + ^RCTdd , (5.7a) 

^bb = ^LCTaa " ^inC^bb + ^R(^ff , (5.7b) 

C^cc = ^in<^bb + i^io'cd - CFdc) + ^ RC^gg " ^ LC^ cc , (5.7c) 

^dd = -^RC^dd + i^icTdc - CTcd) - ^LCTdd , (5.7d) 

(Tee = TlO-cc + i^{(Tef - afe) + T^CTm , (5.7e) 

^ff = ^LCTdd - ^R(^ff + i^{(^fe - CTef) " ^inC^ff , (5.7f) 

&gg = Tir,aff -TRagg-TLdgg , (5.7g) 

^hh = ^LCTgg - ^Rf^hh , (5.7h) 

&cd = iiE3 - E^)(7,d + i^{(Jcc - (Jdd) - l/2(2ri + VR)(7,d , (5.71) 
(7e/ = i(-E3 - E2)aef + iVt{a^e - CTff) - l/2(rj„ + ri^)(7e/ + Tiacd , (5.7j) 

and the resonant current flowing through this system is I/e = Tuladd + cr// + CTgg + CThh]- 

VI. GENERAL CASE 

Now utihzing the results obtained in the previous sections we can write the rate equations 
for the general case. These equations describing the time evolution of the density matrix 
aab{t) of the device are as follows: 

^aa^i 51 ^abicTab " CTfea) " CTaa X] ^'^^d. + X] ^cc^c-^a , (6.1a) 
b{^a) d{^a) c(^a) 

^ab — K^b — Ea)(Tab + M X] ^ab'^b'b — ^ ^aa'f^a'b 

-\(^ab I r„^d+ ^b^d\ + \ aa'b' {^a'^a + ^b'^b) , (6.1b) 

\d{^a) d{^b) I a'b'^ab 

where flab denote the couplings between non-orthogonal isolated states, as for instance be- 
tween the levels in adjacent wells, and a^a — (Jab- The width Ta-^b is the probability per unit 
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time for the system to make a transition from the state \a) to the state \b) of the device due 
to the tunnehng to (or from) the reservoirs, or due to interaction with the phonon bath, or 
any other interaction, generated by a continuum state medium. Notice that Eq. ( p.la| ) for 



diagonal elements has a classical rate equation form, except for the first term. This term 
describes transitions between isolated states through the coupling with non-diagonal terms. 
Therefore it is responsible for coherent (quantum) effects in the transport. 

The non-diagonal matrix elements are described by Eq. (|6.1b| ), which resembles the 
corresponding Bloch equation, supplemented with an additional term. The latter appears 
whenever a one- e/ectron transition converts the state \a') into \a) and the state \b') into \b). 
The positive sign of the additional term calls forth a suspicion that Eq. (|6.1b| ) might have 



unbounded solutions. This is not the case, however, since any positive contribution from the 
additional term in the equation for &ab has its negative counterpart in the equation for &a'b', 
which corresponds to the conversion {a'b') (ab) and originates from the third term in the 
rhs of Eq. ( |6.1b| ). Moreover, the negative contributions, corresponding to the conversions 
(ab) (a'b) and (ab) {ab') have no positive counterparts. Therefore the coupling with 
continuum leads to negative total balance, and hence, to the damping of non-diagonal matrix 
elements. 

The current through the mesoscopic device is the time derivative of the total charge 
accumulated in the collector. We find that the current is totally determined through the 
diagonal elements of the density matrix of the device by the following relation 

m = eY.cTUt)T^^\ (6.2) 

c 

where |c) are the occupied states in the well adjacent to the collector, and is the partial 
width of the state |c) due to tunneling to the collector. 

Although the non-diagonal density-matrix elements do not enter explicitly in Eq. ( |6.2D , 
they are coupled with diagonal matrix elements in the rate equations (the first term in 
( |6.1a| )), and therefore influence the resonant current. The coupling with non-diagonal ele- 
ments always appears in the rate equation, whenever a carrier jumps from one to another 
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isolated states inside the device. In the absence of such transition as, for instance, in res- 
onant tunnehng through a single well, the diagonal and non-diagonal matrix elements are 
decoupled and the evolution of diagonal density-matrix elements is described by the classical 
rate equation. 

Hence, the distinction between isolated and continuum states becomes very essential in 
the description of quantum transport. At first sight, it may seem that in a real situation 
such a distinction can hardly be carried out, since there are no pure isolated states. For 
instance, a single electron state inside the device is always coupled with the continuum 
states of phonons. However, the corresponding density of states would display peaks in 
energy dependence, and they can be considered as isolated states. Indeed, if we have written 
equations like Eqs. (|2.5| ), (|3.3| ), etc. for such a system, the contribution from these peaks 
in the integrals over continuum states would generate a coupling with non-diagonal density 
matrix elements in the rate equations, just as in a transition between two isolated states, 
Eqs. (H). 

As an example of application of Eq. ( |6.1| ) in the case of strong Coulomb blockade, we 
consider the system shown in Fig. 5. The wells may represent two coupled dots. An electron 
tunnels from the emitter to the first well, and then to the second well into the upper level 
E2. After that it can either relax inelastically into the lower level due to interaction with 
the phonon bath, and then tunnel into the collector, or tunnel out directly into the collector. 
Here Fj^ and F'^ are the partial widths of the upper level, E2, due to coupling to phonon 
reservoir and the collector, respectively, and F^ is the width of the level E^ due to coupling 
to the collector. Let's assume that the Coulomb blockade prevents the system from being 
occupied by two electrons, even in different wells. Then there are four possible states of the 
device \a) - all the levels -^1,2,3 are empty; \b) - the level Ei is occupied; |c) - the level E2 is 
occupied; \d) - the level E^ is occupied. It is clear that the density matrix elements for an 
electron with spin up and spin down inside the system are equal, abb-] = crbbi = (Jbh, and the 
same holds for cXcc and add- Hence, Eqs. ( |6.1| ) can be written in this case as 
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Gaa = -2VL(Taa + 2r'^(Tcc + '2'TRadd , (6.3a) 

dhb = i^o{o-bc - 0-cb) + '^LO-aa , (6.3b) 

&CC = -i^o{o-bc - cTcb) - (r'^ + Tin)acc , (6.3c) 

O'dd = -^R(^dd + ^inO'cc , (6.3d) 

cTbc = i{E2 - Ei)abc + i^oicTbb - cFcc) - ^(r^ + Tin)(Tbc , (6.3e) 



and the dc current /, Eq. (|6.2|) , is given by 

I/e = 2accT'R + 2addTR. (6.4) 

Using aaa + 2abb + 2a cc + 2a dd = 1, Eqs. ( |6.3|) can be easily solved for t ^ oo, yielding for 
the dc current 



(r.n + ry2 2 (r.„ + ry2 

This result shows very peculiar dependence of the dc current of the inelastic width Fj^. 
One could expect, at least for T'j^ ~ Tr, that the current should increase when Tin grows. 
However, as follows from Eq. ( |6.3|) , the current / ^ for Tin — > oo (cf. with another example 
in [|l3). In fact, such an unexpected behavior of the dc current would always take place in the 
presence of coherent transitions between isolated states in carrier transport. For instance, 
it can be traced even in a more simple case of the resonant tunneling through a double well 
structure, Eq. (|4.19| ). One finds that / — when F^ oo. This phenomenon can be 
understood by analyzing Eq. ( |6.1b|) for non-diagonal density matrix elements. In contrast 



with the rate equation for diagonal matrix elements, Eq. ( |6.1a|) , the coupling with continuum 
states always leads to damping of non-diagonal matrix elements. Since the transport through 
isolated states goes only via non-diagonal density- matrix elements, Eq. ( |6.1a| ), the total 



current would always decrease with the growth of the corresponding partial widths. 
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VII. SUMMARY 



In this paper we have studied quantum transport in mesoscopic systems (quantum dots) 
containing finite number of isolated quantum states. Starting with the many-particle wave 
function in the occupation number representation, and integrating out the continuum states, 
we have found the equations of motion for the density submatrix of the system. These equa- 
tions have a form of the master (rate) equations for diagonal density matrix element. But in 
addition, non-diagonal density matrix elements, responsible for transitions between isolated 
quantum states, appear in these equations. If, however, these transitions are generated by a 
continuum states medium, the diagonal and non-diagonal density matrix elements become 
decoupled, and the quantum transport is described by classical rate equations. 

It follows from our derivation that the reduction of many-body Schrodinger equation to 
the modified rate equations for density submatrix of the device can be performed only if two 
conditions are met: first, the energy states of the system which carry the resonant transport 
must be inside the bias, Ep — Ep; second, the width of these states is much smaller than 
the bias. If the second condition is satisfied, but the resonant levels of the device are close 
to band edges, our rate equations cannot be derived. Yet, the method still can be used for 
dc current. However, when the bias is less than the level width, the continuum states of 
the reservoir cannot be integrated out in the manner of Sect. 2, and our method cannot be 
applied. 

We have compared some of our results with ones obtained earlier in the literature. For 
example, for the resonant tunneling through a single dot we obtained the same result as 
Glazman and Matveev [Q. In the case of resonant tunneling through double-dot structure 
we found simple analytical expression for dc current under condition of strong Coulomb 
repulsion inside the dots, when no more than one electron can occupy the dots. The obtained 
expression is very close to that found in 0Jl^ . 



As an application of our equations we considered a more complicated case of the resonant 
tunneling in a coupled dot system, where the inelastic process takes place in the course of 
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transport. It was found that the resonant current decreases with the growth of the inelastic 
width. We found that this anomalous behavior always emerges whenever coherent transitions 
are accompanied by inelastic processes. 
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FIGURES 



FIG. 1. Resonant transport through a single quantum well structure. 
FIG. 2. Resonant transport through a double-well structure. 

FIG. 3. Resonant transport through a double-well structure in the presence of inelastic process. 

FIG. 4. All possible electron states of the device, shown in Fig. 3: |a) - all the levels -Ei,2,3 are 
empty; \b) - the upper level, Ei, is occupied; |c) - the lower level, E2, is occupied; \d) - the level 
E3 is occupied; |e) - the levels Ei and E2 are occupied; |/) - the levels Ei and E3 are occupied; 
\g) - the levels E2 and Es are occupied; \h) - all the levels Ei^2,3 are occupied. 

FIG. 5. Resonant transport through a double-well structure in the presence of inelastic process 
with strong Coulomb blockade effects. 
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